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ABSTRACT 

Two populations of minor bodies in the outer Solar System remain particularly elusive: Scattered 
Disk objects and Sedna-like objects. These populations are important dynamical tracers, and un- 
derstanding the details of their spatial- and size-distributions will enhance our understanding of the 
formation and on-going evolution of the Solar System. By using newly-derived limits on the maximum 
heliocentric distances that recent pencil-beam surveys for Trans-Neptunian Objects were sensitive to, 
we determine new upper limits on the total numbers of distant SDOs and Sedna-like objects. While 
generally consistent with populations estimated from wide-area surveys, we show that for magnitude- 
distribution slopes of a > 0.7 — 1.0, these pencil-beam surveys provide stronger upper limits than 
current estimates in literature. 

Subject headings: Kuiper Belt, Trans- neptunian objects 



1. INTRODUCTION 

A number of deep, narrow-angle "pencil-beam" sur- 
veys for distant Solar System objects have been under- 
taken over the past two decades, first appearing in the 
literature in Tyson et al. (1992). These surveys avoid 
the issue of trailing losses in long exposures by taking 
a large number of shorter exposures, predicting the sky 
rate of motion of sources of interest, then compensat- 
ing for this motion in software before finally stacking the 
images. Due to the large number of rates at which im- 
ages must be combined in order to properly compensate 
for the range of motions real objects can have, the or- 
bital range over which this method is applied is limited 
to maintain computational feasibility. 

Parker & Kavelaars (2010 submitted; hereafter P10) 
re-characterize the orbital limits of several published 
pencil-beam surveys and show that these orbital limits 
are poorly characterized in literature. As the re-derived 
maximum heliocentric distances these surveys were sen- 
sitive to ranges from 150 — 400 AU, we find that these sur- 
veys were sensitive to several dynamically interesting yet 
currently poorly constrained populations; namely, Scat- 
tered Disk Objects (SDOs, eg., Trujillo et al. 2000; here- 
after TOO) and Sedna-like objects (SLOs, eg., Schwamb et 
al., 2009; hereafter S09). SDOs are highly-eccentric, non- 
resonant objects that have perihelia that interact with 
Neptune, and may be the source of Jupiter-family comets 
(eg., Duncan & Levison 1997). SLOs are long-period ob- 
jects that have perihelia high enough (> 70 AU) that 
they are out of reach of the giant planets but aphclia 
not distant enough for galactic tides to have significant 
effect on their orbits. Their current orbits require em- 
placement mechanisms that may be linked to dynamics 
of the Solar birth cluster (eg., Brasser et al. 2006), in- 
teractions with "rogue" planets (eg., Gladman & Chan 
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2006), or close stellar passages (eg., Kenyon & Bromley 
2004). Few objects (one, in the case of Sedna) of cither 
population have been discovered, and our understanding 
of these populations is severely hampered by these lim- 
ited samples. Here we use the newly-derived heliocentric 
distance limits from P10 in conjunction with a simple 
survey simulator in order to determine new limits on the 
SDO and SLO populations. 

2. CHARACTERIZATION OF OUTER LIMITS OF 
PREVIOUS SURVEY GRIDS 

In order to characterize the limits of parameter space 
searched by a survey, P10 generated a large sample 
of synthetic orbits that fall on the imaged field dur- 
ing the time of observations. These synthetic orbits 
spanned large ranges of heliocentric distance (20 — 500 
AU), inclination (0° - 180°), and eccentricity (0-0.999). 
Ephemeris software was used to estimate the on-sky mo- 
tion of each synthetic orbit generated, and then test the 
on-sky motion to ensure that it is within the ranges 
searched in the original survey. This characterization 
was performed for those surveys that included enough 
information in the literature to accurately reproduce the 
original search-grids. 

Three surveys were characterized: Fraser et al. (2008), 
Fraser & Kavelaars (2009), and Fuentes et al. (2009). 
Each survey was found to be sensitive to distances con- 
siderably greater than claimed, and this limit varied 
somewhat with inclination. Additionally, each survey 
was sensitive to inclinations as high as 180°, even though 
the target maximum inclinations were 70° — 90°. Table 1 
contains the derived limits and properties of each survey. 

While P10 did not characterize Bernstein et al. (2004; 
hereafter B04), we include it for analysis here. The au- 
thors claim that their survey is sensitive to motions of 
sources at distances greater than 1000 AU and to incli- 
nations as high as 45°. Because B04 imaged fields near 
quadrature (while the three surveys above were taken 
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Table 1 

Pencil-Beam Survey Characterization 

dmax (i = 0°) d max (i = 180°) Am ~ 

Survey (AU) (AU) (mag) (deg 2 ) 

Fraser et al. 2008 164 184 25.25 a 3 

Fraser & Kavelaars 2009 360 390 26.76 0.255 

Fuentcs ct al. 2009 220 245 26.86 0.255 

Bernstein ct al. 2004 1000 b NA 28.5 0.019 
a : Area- weighted mean R$q. 
b : Not characterized by Parker & Kavelaars 2010. 



near opposition), it is not straightforward to apply the 
results of P10 to this survey. However, we are confident 
that adopting 1000 AU as a conservative distance limit 
for B04 is an adequate first-order estimate. 

3. UPPER LIMITS ON DISTANT POPULATIONS FROM 
PREVIOUS SURVEYS 

3.1. Survey Simulations 

With the characterized outer limits of several litera- 
ture pencil-beam surveys for TNOs, we can explore the 
upper limits these surveys put on distant populations. 
Few detections from these surveys are followed-up after 
enough time has elapsed for orbital properties other than 
heliocentric distance and inclination to be measured with 
any precision. With these short-arc measurements, it is 
difficult to identify the population any individual object 
belongs to unless it can be constrained by d and i alone. 
By determining the most distant detection in each sur- 
vey under consideration and simulating a population that 
spends at least some fraction of time outside of this most 
distant detection, we can determine an upper limit for 
the size of this population based on the non-detection 
of any sources at greater distances. Two such distant 
populations are the SDOs and SLOs. 

We simulate each population by sampling objects from 
distributions of heliocentric distance d and absolute mag- 
nitude H. For our model populations, we sample H from 
a power-law absolute magnitude cumulative-distribution 
function (CDF) with the form 

N(< H) = N< Hmax W a{H ~ Hma °' ) , (1) 

where -/V<# max is scaled to produce the total number 
of objects in our sample, given H max is the maximum 
absolute magnitude of interest. 

Because the magnitude distributions of TNOs are 
thought to break to a lower slope at a characteristic 
magnitude and therefore depart from a true power-law, 
we will model our populations with magnitudes brighter 
than H ~ 10 to avoid the putative break. Since we do 
not know the slope of the magnitude distributions a pri- 
ori, we will test our models over a plausible range of a. 

For SDOs, heliocentric distance is sampled from two 
distributions. The first is a differential radial distribution 
with the form 

d/Voc oT 1 - 5 ^ (2) 

with the exponent of —1.5 as used in TOO when per- 
forming similar estimates of distant populations. In the 
TOO models, synthetic SDOs were generated over dis- 
tances of 34 — 366 AU. For our estimate of the upper 
limit of this population, we use the same distance range. 
The second distribution is uniform in heliocentric dis- 
tance over the same distance range, in line with the uni- 



form distribution measured by Kavelaars et al. (2008) 
and returned by the dynamical simulations of Volk & 
Malhotra (2008). 

We also model the SLO population with two radial 
distributions. The first is the same radial distribution 
defined in Eqn. [2 but over 70 - 1000 AU. The second 
distribution is similar to that used by S09, which is the 
intrinsic radial distribution of an orbit with Sedna's a 
and e. The radial distributions used by our simulations 
are illustrated in Figure [TJ with the outer distance limits 
of the surveys under consideration marked. 

We generate a large number of objects drawn simul- 
taneously from each distribution and for each object de- 
termine its apparent magnitude m ~ H + 10 \og(d). For 
each survey, we then determine the following: 

1. Is the object within the heliocentric distance limit 
of the survey? 

2. Given the object's apparent magnitude m, a uni- 
form random probability p, and the survey's effi- 
ciency function rj(m) (clipped at 15%), is the object 
detected (p < 77(771))? 

If an object satisfies both conditions, it is counted as 
"detected" by that survey. We also note whether or not 
the detection is outside the most distant H < 10 detec- 
tion made by the survey. We then scale the number of 
detections in each survey by the fractional area fa of 
the entire sky the survey searched. Given this scaled de- 
tection count for each survey, we determine the mutual 
probability of non-detection of the synthetic population 
more distant than the most distant real detection for all 
surveys. We then determine the initial number of objects 
in the population required to return a mutual probabil- 
ity of non-detection 0.27% (3cr upper limit), given this 
population model. For every a we simulate, we take 
the corresponding derived upper limit on the population 
(jV(< H = 10)) and convert it to N(D < 100 km) by 
assuming a 4% albedo. 

By making detection a function only of H and d, we are 
effectively modeling a population that is isotropic on the 
sky, which implies an inclination distribution of p(i) cx 
sin(i). For a population with an unknown inclination 
distribution this is as reasonable an assumption as any, 
and it significantly simplifies analysis. We note that TOO 
only considered inclinations as high as 35° for SDOs, 
and S09 modeled SLOs with an inclination probability 
distribution that was effectively for i > 35°. In order 
to compare to the earlier estimates on the upper limits 
of the SDO and SLO populations, we scale fn in our 
simulations by a factor of 1/ sin(35°), which has the effect 
of limiting the total area of the sky considered to that 
which inclinations < 35° can populate. 
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Fig. 1. — Cumulative radial distributions and limits used in our simulations. Grey lines: d 1,5 (solid line) uniform (dashed line) radial 
distributions used for SDO models. Black lines: d~ l j distribution (solid line) for SLOs, and intrinsic radial distribution (dashed line) of 
object with Sedna's a and e. Vertical bars: Distance limits for characterized surveys (widths represent variation of limit with i). Vertical 
dotted line: Estimated distance limit for B04. 



Since the field locations of the pencil-beam surveys 
considered here all have ecliptic latitudes near ~ 0°, the 
upper limits we derive only truly reflect an extrapola- 
tion of the population density on the ecliptic. In order 
to determine the total population limit N' implied by 
our surveys given an arbitrary latitude distribution, one 
may simply multiply our upper limit No by the integral 
of the desired latitude distribution P'(l) (where P'(Q°) 
is normalized to 1) between I of to ir/2 radians, and 
then divide by the integral of our assumed distribution 
(Po(0 = cos(Z) between and 0.611 radians, elsewhere) 
over the same range: 



N' 



N 
0.57 



7T/2 



P'(l)dl. 



For example, assuming a latitude distribution for the 
SDOs similar to that determined by Brown (2001; ap- 
proximated here as a gaussian with a = 0.31 radians) 
results in a correction of our upper limits by a factor of 
~ 0.6 -0.7. 

4. SCATTERED DISK POPULATION LIMITS 

Since the perihelia of the SDO population extend well 
inside the outermost detections of the surveys considered 
here, the non-detection of distant sources is not necessar- 
ily a stronger constraint on the population than the total 
number of detections consistent with being SDOs in the 
survey. To account for this, we also count the number 
{N S i m ) of H < 10 simulated detections returned by each 
survey where the distance at detection of the simulated 
object was less than the distance of any real detection 



by that survey. We assume (as a limiting case) that all 
N rea i detections with H < 10 and d > 34 AU made by 
each survey were SDOtQ. We then determine the proba- 
bility of detecting N rea i given N S i m , and derive an upper 
limit when N s i m is ruled out at the 3cr level. If this up- 
per limit is stronger than the upper limit derived by the 
non-detection of distant objects for a given a, we adopt 
it as our upper limit for that slope. 

In addition to modeling the pencil-beam surveys, we 
have derived separate limits on the SDO population by 
passing the detections of TOO through our simulation. 
TOO had several verified real detections, which allow us 
to derive both upper and lower limits for the population. 
Our upper limits in the d~ 15 radial distribution case 
match those estimated by TOO well, while our lower limits 
are a factor of ~ 2 higher than the values quoted in 
TOO. This discrepancy is likely due to variations between 
the uniform inclination distribution used by TOO and the 
uniform latitude distribution used here. 

Figure [2] illustrates the results of our simulations for 
the i < 35° SDO population, the SDO limits extrapo- 
lated from our simulations of the detections in TOO, and 
limits on the Main Kuiper Belt (MKB) population from 
the "nominal" model of Kavelaars et al. (2009; hereafter 
K09). For increasing a, the pencil-beam survey upper 
limit on the population decreases, while the upper limit 
from TOO increases. For both radial distribution models, 
the point where the two limits cross defines an absolute 
upper limit on the population for any CDF power-law 

1 Because these surveys do not provide accurate measurement of 
orbital parameters other than heliocentric distance and inclination, 
no more refined characterization of individual objects is possible. 
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Fig. 2. — Limits on total population of i < 35° Scattered-Disk Objects with diameters < 100 km vs. CDF slope a. Solid lines: 3cr upper 
limit from pencil-beam surveys considered here (black) and extrapolated upper- and lower-limits from the detections in TOO (light gray) 
using the d~ 1,5 radial distribution. The inflection point in the black line occurs at the value of a where the upper limit constraint on 
the population of SDOs transitions from being driven by detections at small distance (a > 0.5) to being driven by the lack of detections 
at large distance (a < 0.5). Dashed lines: 3cr upper limit from pencil-beam surveys considered here (black) and extrapolated upper- and 
lower-limits from the detections in TOO (light gray) using a uniform radial distribution. Points: limits estimated by TOO, with 3<r limits 
marked by error bars. Gray band: 3cr limits on Main Classical Kuiper Belt population by "nominal" model of K09, with cimkb = 0.72. 



slope a (given that radial distribution model). For the 
d~ 15 radial distribution model, this crossover occurs at 
approximately a = 0.95, with a maximum population 
size of N(D > 100 km) < 3.5 x 10 5 objects. For the uni- 
form distribution model, this crossover occurs at the shal- 
lower slope of approximately a = 0.8, with a maximum 
population size limited to N(D > 100 km) < 2.5 x 10 6 
objects. These limits translate into 5 — 40 times the up- 
per limit population of the K09 "nominal" model of the 
MKB. 

5. SEDNA-LIKE POPULATION LIMITS 

The most distant detections (~ 43 — 65 AU) in each of 
the pencil-surveys considered here are all interior to the 
modeled perihelia of the SLO populations (70 — 76 AU). 
The limiting magnitudes of these surveys allowed them 
to detect objects as faint as H ~ 10 (B04 limit at 70 
AU). While these surveys searched over 3,000 times less 
area than S09, their flux limits are up to nearly seven 
magnitudes fainter than the limit of S09 at perihelion, 
permitting a significant improvement on the limits for 
steep CDF slopes. 

Comparing the SLO population estimates of S09 to the 
limits derived in this work is not straightforward. S09 
was sensitive only to the largest objects in the population 
while the surveys analyzed in this work were sensitive to 
much smaller objects. We have modeled the detection of 
Sedna by S09 with the same populations used to simulate 
the pencil-beam surveys (with both radial distribution 
models) and extrapolated the 3er upper-limits the S09 



survey put on the bright population to diameters as small 
as 100 km. 

Figure[3]illustrates the results of our simulations for the 
SLO population (given both radial distribution models), 
the SLO limits we derive from the detection of Sedna by 
S09 extrapolated to the D ~ 100 km regime, and limits 
on the MKB population from the "nominal" model of 
K09. The large difference (factor of ~ 50) between the 
two limits illustrates the sensitivity to assumed radial 
distribution models. 

The crossover where the upper limit from pencil-beam 
surveys becomes a stronger constraint on the SLO pop- 



ulation is approximately a ~ 0.7 for the d 



radial 



distribution and a ~ 0.63 for the Sedna-like radial dis- 
tribution. At the slope where the limit derived here is 
stronger than the extrapolated S09 result for each radial 
distribution model, our upper limit on the SLO popula- 
tion is ~ 10 times larger than the K09 upper limit on the 
MKB for the distribution, and the largest the SLO 
population can be for any CDF slope with this radial 
distribution is roughly 40 times the MKB upper limit. 
Given the Sedna-like radial distributions, the SLO pop- 
ulation can be up to ^1300 times more populous than 
the upper limit on the MKB (with the peak occurring 
at a ~ 1) before the lack of detections in current pencil- 
beam surveys would be incongruous with expectations. 

5.1. Requirements for Future SLO Surveys 

If a survey is undertaken to detect the SLO population, 
or at least drive down the upper limit on the population, 
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Fig. 3. — Limits on total i < 35° Sedna-like population with diameters < 100 km vs. CDF slope a. Solid lines: 3tr upper limit from 
pencil-beam surveys considered here (black) and extrapolated from S09 (light gray) using the d~ 15 radial distribution. Dashed lines: 
3<t upper limit from pencil-beam surveys considered here (black) and extrapolated from S09 (light gray) using Sedna's intrinsic radial 
distribution. Gray band: 3<r limits on Main Classical Kuiper Belt population by "nominal" model of K09, with cimkb = 0.72. 
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Fig. 4. — Required survey area f2 (square degrees) vs. survey 
R limiting magnitude in order to obtain stronger 3cr upper limits 
on Sedna population than derived here, for magnitude-distribution 
power-law slopes a = 0.6, 0.8, and 1.0. Diamond: expected depth 
and area of a single LSST "deep drilling" field. 

how deep and over what area must the survey search? 
Given the uncertainty in the properties of the population, 
we have estimated the area Q a survey must search, given 
a limiting magnitude R\[ m (and a = 0.6, 0.8, and 1.0), in 



order to obtain stronger 3c upper limits on the SLO pop- 
ulation than those derived here. Figure [4] illustrates the 
results for survey limiting magnitudes of R\[ m = 25 — 28. 
The model survey is assumed to be sensitive to motions 
out to 1000 AU. A single Large Synoptic Survey Tele- 
scope (LSST) "deep drilling" field (Chesley et al. 2009) 
will probe to nearly the same depth as Bernstein et al. 
(2004) but over 500 times the area. If the methods used 
to search this component of the LSST data for TNOs 
is sensitized to distant objects, the chances for finally 
obtaining a statistically useful sample of the Sedna-like 
population are very good. 

6. SUMMARY 

We have used the re-characterization by Parker & 
Kavelaars (2010) of the orbital sensitivity of existing 
pencil-beam surveys in literature to derive new upper 
limits on distant populations of minor bodies in the So- 
lar System. P10 considered three surveys (Fraser et al. 
2008; Fraser & Kavelaars 2009; and Fuentes et al. 2009) 
and found that for each survey the outer limit was sig- 
nificantly more distant than claimed in the original pub- 
lication. Using these newly-derived outer limits, we de- 
termined the constraints these surveys (plus Bernstein 
et al. 2004) put on the Scattered Disk and Sedna-like 
populations. 

The combination of survey types allows us to limit the 
maximum possible size of the population for any CDF 
slope. Our limits constrain the Scattered Disk to be 
less than 5 — 40 times more populous than the Main 
Kuiper Belt, depending on the radial distribution of ma- 
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terial. For most plausible CDF slopes, SDOs are better 
constrained by relatively shallow, wide area surveys like 
those presented in Trujillo et al. (2000). For both a e?~ 15 
and a uniform radial distribution model, we estimate the 
maximum population size at any CDF slope to less than 
N(D > 100 km) < 3.5 - 25 x 10 5 , respectively. 

A strong sensitivity to the assumed radial distribution 
is demonstrated by the Sedna-like population, as up to 
1300 times the D < 100 km Main Kuiper Belt popula- 
tion could reside on Sedna-like orbits (as assumed by the 
limits presented in S09) and still remain undetected in 
existing surveys, whereas if the material follows a d -1 ' 5 
radial distribution this upper limit is reduced to 40 times 
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the current population of the Main Kuiper Belt. 

At present the upper-limits on the population size de- 
termined by lack of detection in existing pencil-beam sur- 
veys do not appear to conflict with the population esti- 
mates based on wide-area, shallow surveys. This implies 
that the lack of distant detections in current-generation 
surveys are not indicative of any incongruences in survey 
sensitivity calibrations or other issues. 
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